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Abstract 

Using fully relativistic GEANT4 simulation tool kit, the transport of ener- 
getic electrons generated in solar flares was Monte-Carlo simulated, and resultant 
bremsstrahlung gamma-ray spectra were calculated. The solar atmosphere was ap- 
proximated by 10 vertically-stacked zones. The simulation took into account two 
important physical processes,that the bremsstrahlung photons emitted by precipitat- 
ing relativistic electrons are strongly forward beamed toward the photosphere, and 
that the majority of these gamma-rays must be Compton back-scattered by the solar 
atmosphere in order to reach the observer. Then, the Compton degradation was found 
to make the observable gamma-ray spectra much softer than is predicted by simple 
analytic calculations. The gamma-ray signals were found to be enhanced by several 
conditions, including a broad pitch-angle distribution of the electrons, a near-limb 
flare longitude, and a significant tilt in the magnetic field lines if the flare longitude is 
rather small. These results successfully explain several important flare properties ob- 
served in the hard X-ray to gamma-ray range, including in particular those obtained 
with Yohkoh. A comparison of the Yohkoh spectrum from a GOES X3.7 class limb 
flare on 1998 November 22, with a simulation assuming a broad electron pitch- angle 
distribution, suggests that gamma-rays from this particular solar flare were a mixture 
of direct bremsstrahlung photons and their Comptonization. 

Key words: solar flare: gamma-ray: Compton scattering 



1 



1. Introduction 



Solar flares are sudden, brief, and powerful outbursts, releasing an energy of 10 28-34 
ergs on time scales of a second to several tens of minutes. In solar flares, both protons and 
electrons are accelerated to non-thermal energies. Although these energetic charged particles 
are strongly affected by the solar magnetic fields and cannot usually be detected, we can utilize 
neutral secondary particles, neutrons and photons, to probe them. Among them, hard X-ray 
to gamma-ray continua, produced via bremsstrahlung process, provide the best diagnostics of 
the flare-accelerated electrons. 

Although high-spatial-resolution images of flares in the hard X-ray and gamma-ray en- 
ergies have become available by the Yohkoh HXT (Kosugi et al. 1991) and RHESSI (Lin et al. 
2004) respectively, we still depend heavily on the spectral data when attempting to study the 
non-thermal electrons. Actually, to calculate back the electron spectrum from the observed 
gamma-ray spectra has been an important issue in solar physics. As we understand, flare- 
accelerated electrons having a power-law distribution penetrate the solar chromosphere, and 
produce bremsstrahlung photons via "thick-target" emission. The thick-target emission is a 
process wherein non-thermal electrons impinge on a thick matter, and lose almost all their 
energies through repeated Coulomb collisions while emitting X-ray/gamma-ray photons via 
bremsstrahlung. 

In the thick-target condition, the spectrum of X-ray emitting electrons is formed by an 
equilibrium between the injection of newly accelerated electrons and the loss of their energies 
predominantly through collisions. If the injected electron number spectrum is given by a power- 
law as F(E) ~ E~ 5 (electrons keV" 1 s" 1 ), with E the electron kinetic energy and 5 a constant 
called electron index, then we expect the X-ray photon-number spectrum also to take a power- 
law from as 



where 7 is a quantity called photon index. In 1974, Brown found an approximate relation as 



in non-relativistic regime, assuming that all emitted photons are collected. The photon spec- 
trum thus becomes relatively flat, because of "loss-flattening" effect working on the electrons; 
as they penetrate deeper into the target material, the spectrum gradually hardens because 
lower-energy electrons have shorter time scales of energy loss through Coulomb collisions. 

Besides the thick-target condition, another extreme case called thin-target condition 
(Brown 1971) was studied extensively. This represents a situation wherein a bunch of energetic 
electrons pass through a thin medium while emitting bremsstrahlung photons, but they leave 
the emission region with their energies almost unchanged. In the thin-target condition, the 
energy distribution of the incident electron is not significantly affected by collisions with the 




(1) 



7 = 5-1 



(2) 



2 



target. The emergent photon spectrum is steeper, because the electrons are free from the loss 
flattening effect and is characterized as (Brown 1971) 

7 = 5 + 1/2. (3) 

Since then, there have been many attempts to improve equations (2) and (3), particularly 
toward relativistic regime. Some utilized purely analytic calculations of various probability dis- 
tributions, while others employed numerical solutions to the Fokker-Planck equation, or Monte 
Carlo techniques. In electron energies up to 10 MeV, Bai (1982) investigated angular depen- 
dence of the bremsstrahlung in the hard X-ray range, and discussed the associated electron 
transport. Murphy et al. (1987) dealt with the particle transport by employing a simple thick 
target emission model, assuming an isotropic electron distribution. Dermer & Ramaty (1986) 
conducted a detailed study of the directional bremsstrahlung emission by anisotropically accel- 
erated electrons. Miller & Ramaty (1989) made a calculation including pitch angle scattering 
of electrons by plasma waves. These previous works have indicated that the gamma-ray con- 
tinua of solar flares can be generally explained in terms of chromospheric thick-target emission, 
produced by electrons having a power-law energy distribution with 2 < 5 < 4. 

In spite of these extensive studies, the estimation of electron spectra from flare gamma- 
ray data would not be self-contained unless we take into account another important physical 
process; namely, Compton scattering. The relativistic bremsstrahlung from precipitating elec- 
trons must be significantly forward-beamed, so that a majority of gamma-ray photons must be 
Compton back-scattered by the solar materials in order to reach us. This process will strongly 
degrade the gamma-ray energies, and affect the photon spectrum to be observed. To treat this 
kind of physical condition which involve multiple Compton scatterings, we clearly need to use 
Monte Carlo techniques. 

The effect of Compton scattering has been studied by several authors, beginning, e.g., 
with Bai & Ramaty (1978) who used Monte-Carlo simulations in such early days. Kontar 
et al. (2006) considered this issue using a semi-analytic Green's function by Magdzias & 
Zdziarski (1995), and calculated "primary" photon spectra from observed ones. These works 
have clarified that the Compton back-scattering can significantly modify the bremsstrahlung 
spectra. However, these pioneering works have treated electrons with sub- MeV kinetic energies, 
leaving higher energies less studied. 

In the present paper, we calculate not only the electron transport and bremsstrahlung 
emission, but also the photon propagation via Compton process. For this purpose, we employ 
a Monte Carlo simulation toolkit named Geant4, which is widely used in experimental high 
energy physics. Our results imply that the Compton scattering has significant effects on the 
observed gamma-ray spectra. Together with full- relativistic effects, this is expected to require 
revision to the relation such as equation(2) and equation(3). We show the relevance of these 
results to actual flares, by comparing the simulated gamma-ray spectra with those observed 
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with Yohkoh from a sample of solar flares. 
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2. Method of Simulation 



2.1. The Geant4 toolkit 

The Geant4 Monte Carlo simulation toolkit (Agostinelli et al. 2003), developed in ex- 
perimental high energy physicists, has four important features. Firstly, it allows a very flexible 
construction of "geometry" in which particles interact. Secondly, it can "track" individual par- 
ticles while monitoring their physical quantities such as the energy, momentum, and position; 
the name Geant in fact comes from "GEometry ANd Tracking". Thirdly, it allows us to in- 
corporate the desired interactions; in the present case, Coulomb scattering and bremsstrahlung 
for electrons, and Compton scattering for photons. Finally, the toolkit has been tested and 
calibrated extensively by many investigators under different conditions. 

2.2. Particle injection, interaction, and tracking 

In our Monte Carlo simulation using Geant4, we basically shoot a parallel beam of 
a large number of electrons, either vertically or slantly, into our "solar atmosphere" to be 
explained later. The electrons are assumed either to be mono-energetic, or to have a power-law 
energy distribution with an index 8. In the latter case, the initial electron energies are assumed 
to be in the range of 1-100 MeV, so as to simulate gamma-ray flares. 

We implemented the following physics processes which we require our particles to obey in 
our simulation: ionization, multiple scattering, and bremsstrahlung for electrons, while photo- 
absorption, Compton scattering and electron-positron pair creation for photons. We suppressed 
the production of secondary electrons for simplicity. Because those secondary electrons (typi- 
cally less than 100 keV) have much lower energies than the primary ones, they can be considered 
negligible. 

One of the essential features of Geant4 is step-by-step tracking of each individual 
particle, considering all the implemented elementary processes. The particle tracking algorithm 
roughly consists of the following 4 steps. 1) The initial particle velocity is calculated. 2) 
Each physical process calculates the mean-free-path of the specified particle under the relevant 
interaction, and generates a random number around the mean. Then it proposes this random 
number as a step length. 3) Among various competing processes, the interaction which proposes 
the shortest physical length is adopted. 4) Track properties, such as the kinetic energy and 
momentum, of the current particle are updated. 

Among the implemented processes, the Coulomb scattering of electrons need special 
consideration, because its mean free path is much smaller than those of the other processes 
(e.g., ~ 10~ 2 of that of bremsstrahlung for a 10 MeV electron). As a result, its full Monte Carlo 
treatment would demand too much computing times. Therefore, the expected total energy 
loss of an electron due to multiple Coulomb scatterings, in each step which is determined by 
the processes other than the Coulomb scattering itself, is computed based on semi-empirical 
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formulae by Messel & Crawford (1970). Fluctuations around the mean energy loss is expressed 
by generating random numbers after Lassila-Perini & Uarban (1995). Similarly, the total spatial 
displacement caused by the multiple scattering in each step is calculated based on Lewis method 
(Lewis 1950), which solves a diffusion equation of electrons in the matter. Fluctuations around 
the expected mean is again represented by random numbers (Agostinelli et al. 2003). Thus, 
Geant4 simulates the effect of Coulomb interactions after a given step as a statistical sum of 
a large number of scatterings, instead of treating them one by one. 

In this way, the injected electrons and the bremsstrahlung-generated photons are tracked, 
until they escape out of the model boundary, or they become less energetic than a certain 
threshold energy which we set at 1 keV. We collect these photons under a specified condition, 
and make their spectra. 

The simulation of flare-accelerated electrons, precipitating onto the chromosphere, re- 
quires yet another physical process to be implemented: electron gyration around the solar 
magnetic fields. Its exact treatment, however, would make the simulation extremely time- 
consuming, because the gyration radius of a ~ 100 MeV electron in a ~ 100 Gauss field is as 
small as ~ 10 m. Therefore, we have decided not to take gyration into our simulation, but 
instead, to consider slant injections; as illustrated in figure 1, injecting electrons at a polar 
angle a, and collecting the emergent photons at another polar angle /3 but all over the azimuth 
angle tp, are equivalent to simulating electrons gyrating with a constant pitch angle a, as long 
as the field lines are close to normal to the photosphere. The case of tilted magnetic fields can 
be reproduced by a superposition of different injection angles, a. 

2.3. Electron transport under a simple geometry 

Before actually simulating the electron transport and gamma-ray emission employing a 
realistic model for the solar atmosphere, we carried out simpler simulations with two purposes 
in mind; to validate our simulation code by examining whether each of the relevant physical 
processes is correctly reproduced, and to grasp the essence of the physics to be investigated. We 
accordingly approximated the sun as a hydrogen gas box of 10,000 km x 20,000 km x 20,000 km, 
having a uniform density of 3.2xl0~ 7 g cm -2 which is the density at the solar photosphere. 
Hereafter we call this model "simple solar atmosphere model" ; it is one of three uniform target 
models employed in the present paper. 

We vertically injected 10 5 mono-energetic electrons with an initial kinetic energy of 50 
MeV into our "simple solar atmosphere" . In order to monitor how the Coulomb interaction af- 
fects the spectrum and angular distribution of the electrons, we prepared imaginary boundaries 
at different depths, 100, 200, and 300 km from the injection surface, and collected electrons 
at each boundary. Figure 2 a shows the electron spectra obtained in this way, while Figure 
2b the associated angular distributions. From these, we can see that the maximum energy of 
the spectrum decreases as the electrons penetrate deeper. Since a 100 km thick slab in our 
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"solar atmosphere" has a column density of 3.2 g cm -2 , an electron of energy ~ 50 MeV pass- 
ing through it is expected to lose ~ 15 MeV, as we can easily calculate from the Bethe-Bloch 
formula. The results in Figure 2a indeed meet this expectation. The geometrical thickness 
(10,000 km) of the model is thus large enough to make it physically thick to Coulomb loss, for 
electrons with initial energies up to ~1.5 GeV. Furthermore, the electron energy distribution 
is seen to broaden as the energy loss proceeds. 

Figure 2b reveals that the injected electrons, initially forming a parallel beam, are grad- 
ually deflected due to Coulomb interactions with ambient atoms. The simulated angular distri- 
butions approximately consist of two components: a Gaussian component caused by small-angle 
multiple scatterings, and a power-law tail due to large-angle Rutherford scatterings. 

In a similar way, we simulated the transport of energetic electrons which has a power-law 
distribution over the 1-100 MeV range with 5 = 1.2. The spectra and angular distributions, 
measured at depths of 100, 200, 300, 400, and 500 km, are presented in figure 3. The maximum 
electron energy (initially at 100 MeV) thus decreases according to the Bethe-Bloch formula. 
Furthermore, the low-energy part of the spectrum gradually flattens, because the energy loss 
per unit length is roughly independent of the electron energy, and hence lower-energy electrons 
suffer a larger fractional energy loss than higher-energy ones. This is essentially the same as 
the "loss flattening" effect mentioned in section 1. 

The electron index we employed here, 5 = 1.2, is significantly smaller (flatter) than would 
be achieved in the standard diffusive shock acceleration process (Blandford & Ostriker 1978; Bell 
1978). This is motivated by two reasons. One is that the effects of Compton scattering are more 
clearly observed as S gets smaller. The other is that such a flat electron distribution could be 
realized in actual solar flares, via, e.g., direct electric acceleration. Hereafter, we hence consider 
both flat (5 < 2) and steep (5 > 2) electron distributions. 

2-4- Bremsstrahlung spectra 

We investigated photon spectra created by the bremsstrahlung process under the two 
representative conditions, namely the thick-target and thin-target conditions explained in sec- 
tion 1. The GEANT4-simulated gamma-ray spectra are compared with analytic formulae, which 
are fully relativistic unlike the non-relativistic calculations of eq.(2) and eq.(3). 
2.4-1- Thin-target emission 

To validate the bremsstrahlung process under the thin-target condition which is simpler 
of the two, we let 10 s mono-energetic electrons vertically penetrate a second simplified model 
atmosphere, namely a uniform hydrogen gas with a size of (1000 km) 3 , having a density of 
1.0 x 10~ n g cm -3 ; this is henceforth called "thin atmosphere model". The model has a 
column density of 1.0 mg cm -2 , so that the injected electrons lose their energies by no larger 
than 1 %, in agreement with the thin-target condition. We collected the information of each 
bremsstrahlung photon when it was created, and obtained the spectra shown in figure 4. 
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Also shown in figure 4 are analytically calculated spectra, using a fully relativistic formula 
(Schiff 1951; Koch & Motz 1959), as 
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Here, da^jdk is a cross section for an electron of total energy eom e c 2 (including the rest mass 
energy) to interact with a target atom of charge Z, and emit a bremsstrahlung photon of energy 
km e c 2 , thus achieving a final energy of em e c 2 = (e — k)m e c 2 . The classical electron radius is 
denoted r e . 

As seen in figure 4, the difference between our simulation and the analytical calculation 
is less than 10 %, indicating that the bremsstrahlung process has been correctly implemented. 
In the fully relativistic regime(eo 3> l,e ^> 1), in particular, the emitted gamma-ray spectrum has 
a power-law shape with a photon index 7 ~ 1.0, or an energy index of ~ 0, up to the maximum 
photon energy that is identical to the initial photon energy. This is reasonable, because eq.(4) 
predicts da^jdk oc k~ l when e ~ e and hence b^$> 1. 

Figure 5 shows the bremsstrahlung photon spectra calculated under the same thin-target 
condition, but when the electrons are assumed to have a power-law spectrum distributed from 1 
to 100 MeV, with 5 = 0.8 or 1.5. The simulated spectra again exhibit approximately power-law 
like shapes, with a steeper slope than in the case of mono-energetic electrons because softer 
electrons can obviously emit only lower-energy photons. 

In figure 5, we also show photon spectra calculated analytically using equation (4) as 

N(k)=rii v(e )J e (e )—de , (6) 

Jk+l UK 

where N(k) is the photon number spectrum, J e (eo) is the injected electron spectrum (oc eg S in 
the present case), rij is the target ion density, and v(eo) is the velocity of an electron having an 
energy eo- The simulation agrees with the analytic results within ~ 20%. Thus, the Geant4 
simulation and the fully-relativistic analytic formalism consistently indicate 7 ~ 1.3 for 8 — 0.8 
and 7 ~ 1.6 for 5 = 1.5, both measured in a typical gamma-ray energy range of 1 to 10 MeV. 
These photon indicies are significantly smaller (flatter), particularly for 5 = 1.5, when compared 
to the conventional non-relativistic relation of equation (3). 

In this way, we repeated the Geant4 simulation with difference values of the electron 
index, and calculated the Bremsstrahlung photon index over an energy of 1-10 MeV. The 
obtained relation between the electron index and the photon index is shown in Figure 6, in 
comparison with the prediction by equation (3). We see that the formula of equation (3) in the 
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non-relativistic regime, namely 7 = <5 + 1/2, is no longer valid in the present relativistic regime, 
with the discrepancy increasing toward larger electron indicies. 

The cross section of bremsstrahlung as a function of emitted photon energy and angle 
is given by Schiff (1951) and Koch & Motz (1959) as 

el + e 2 4y 2 e 
W + lfel (y 2 + l)V 
where "reduced photon angle" y, together with its function M(y), is given by 

and 8q is the initial angle of an emitted photon. Figure 7 shows equation (7) as a function 
of eo#o, where the results of the simulation using the thin atmosphere model are superposed. 
Thus, the whole emission becomes concentrated roughly within eo#o < 1; the emission is strongly 
forward-peaked to within |6> | < Cq 1 . 
2.4-2. Thick-target emission 

The thick-target emission is significantly more complex than the thin-target case, be- 
cause it results from electrons while their energies change continuously due to Coulomb scatter- 
ing. Moreover, it has a higher practical importance, because the flare hard X-rays and gamma- 
rays are thought to be produced mainly through this process in the solar chromospheric and 
denser regions. 

Under the thick-target condition, we first investigated the emission from mono-energetic 
electrons, as we did for the thin-target case. We shot 10 7 electrons with an initial kinetic 
energy of 10 MeV or 50 MeV, into a hydrogen target of which the size is (10 4 km) 3 ; this is our 
third simplified model atmosphere called "thick atmosphere model" . Figure 8 shows the spectra 
produced by collecting all photons emitted by each electron while its energy gradually decreases 
down to ^ 1 keV. Compared with the thin-target emission from mono-energetic electrons, the 
photon spectra are somewhat steeper, and roll over more prominently toward the maximum 
energy. Obviously, this is because the electrons keep losing their energies while radiating. 

In figure 8, we also show the results of analytic calculations, which are formulated as 

n *le=fe+i da 



rt\t=k+i da 
N(e ini ,k) = Hi / v TT dt 
Jt\ e=H . dk 



f e ^i ( de\ 1 da 

=ni L v {dt) dk de ' (9) 

where t is the time, e- in \ is the initial total energy of an electron, de/dt is its energy loss rate, v 
is the instantaneous electron velocity, rij is the ion density of the target plasma, and % refers 
to equation (4). Again, the simulation and analytic calculation agrees within < 20%. 

We investigated the thick-target emission from power-law distributed electrons, with 
5 = 0.8 or 1.5, over the 1-100 MeV energy range. Employing 10 6 electrons, we collected infor- 
mation of the photons in the same way as in the preceding simulations, and obtained figure 9. 
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Designating the electron injection rate (electrons cm 2 s x ) by /(e^) and the target area by A, 
and employing equation (4), the photon production rate is analytically given by 

/^max 
N(e- mh k)f(e mi )de- n 
-min 

f £ max /-Cini d<J 

= A / rfe ini /(e in i) / de—v(e) 

Je min Jk+1 ak 

This prediction is again show on the figure. The simulation agrees within ~ 20% with the 
analytic calculation via eq.(10). 

As can be seen from figure 9, the photon spectra emitted by power-law distributed 
electrons under the thick-target condition again exhibit power-law shape up to ~ 10 MeV, like 
in the case of the thin-target emission (figure 5). The obtained photon index is 7 ~ 1.4 for 
5 ~ 0.8 and 7 ~ 1.5 for 5 ~ 1.5. In contrast, the hardest portion of the spectra exhibits rather 
opposite trend, and cannot become as hard as implied by equation. (2), since the maximum 
energy of electrons decreases as they propagate deeper through the thick target. 

Figure 6 (right) is the same comparison as figure 6 (left), but for the thick-target con- 
dition. Photon indicies in this case are generally smaller compared with thin target emission, 
because of the loss-flattening effect (section 1). In addition, 7 depends more weakly on 5 than 
in the thin-target case. 



de 
dt 



(10) 
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3. Solar Simulation 



3.1. A more realistic model atmosphere 

In the previous section, we treated the solar atmosphere as a uniform gas slab. This 
is basically reasonable, because the energy losses of electrons and photons are primary deter- 
mined by the column density of matter along its path. The actual solar atmosphere, with the 
strong vertical density gradient, could be transformed into such a uniform zone, as long as all 
the particles are traveling vertically. However, particle trajectories with significant transverse 
components can no longer be conserved by such a transform, and hence we need to construct 
a more realistic solar model taking into account its vertical density gradients. We here refer to 
"Harvard-Smithsonian reference atmosphere" (Gingerich et al. 1971) for the solar mass-density 
profile, and approximate this numerical table analytically by an exponential function as 



Here, z is the vertical coordinate measured from the photosphere, p(z) is the mass density in 
the unit of g cm -3 , and h is the scale height which is ~ 400 km for z < (the solar interior) 
and ~ 110 km for z > (the coronal region). 

We construct our new solar atmosphere model as a series of boxes (or "zones" ) stacked 
in the vertical direction. Individual zones are defined to have different densities but to share 
approximately the same column density, and the density within each of them is assumed to be 
uniform. Each zone is assumed to have a lateral extent of 2 x 10 4 km, to simulate the cross 
section of a typical magnetic loop. On this transverse scale, the curvature of the solar surface 
is negligible. As illustrated in figure 10, we divide the coronal (z > 0) region into 5 zones, each 
having approximately the same column density of 0.70 g cm -2 . 

According to the Harvard Smithonian model, the total column density above z = is 
3.5 g cm -2 , which corresponds to a depth of ~ 100 km in our simple atmosphere model. Also, 
the sum of the 5 coronal zones in our new model gives an overall column density of 3.5 g cm" 2 , 
thus faithfully representing the actual solar atmosphere for z > 0. While traversing each zone, 
a 100 MeV electron loses ~ 15 MeV in total , of which only ~ 20 % is in radiation whle the 
rest is in Coulombic loss. 

In simulating the solar interior {z < 0), we use 10 thicker zones having 3.5 g cm -2 each; 
again, this corresponds to a depth of ~ 100 km in our simple atmosphere model. With the 
overall column density of 35 g cm -2 below z = 0, our new model goes down to a depth of ~ 500 
km, below the solar photosphere in comparison with equation (11). This depth is considered 
sufficient, because a 100 MeV electron, injected at the top of our model, would lose ~ 100 
% of its initial energy by the time it reaches the bottom of the deepest zone. Furthermore, 
the thicker 10 zones, when summed up, has a Compton optical depth of ~ 10, for a vertically 
precipitating photons; few photons would penetrate deeper. 




(11) 
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Although our "solar atmosphere" is comprised of pure hydrogen, the actual solar atmo- 
sphere contains helium to ~ 10% by number, or ~ 25% by mass. Since the ionization loss of 
electrons is proportional to the electron number density, the helium would decrease the ioniza- 
tion loss per unit mass column density by ~ 10%. Bremsstrahlung per unit mass should remain 
the same, because a helium atom has 4 times higher mass than a proton but 4 times higher 
bremsstrahlung cross section at the same time. The Compton scattering depends only on the 
number of electrons, and hence would decrease by ~ 10% if helium was included. 

3.2. vertical injection 

To examine the flare photon spectra, we vertically injected 10 6 electrons into our solar 
atmospheric model from a height of z = 10, 000 km, and tracked them from that height to 
—500 km, as they lose energy via Coulomb scattering and emit bremsstrahlung photons. The 
electrons are assumed to be distributed initially in the 1-100 MeV range with a power-law 
index 5 = 1.2. The electrons, initially moving in the vertical direction, gets gradually deflected 
as indicated in figure 3; in actual configuration, this corresponds to pitch-angle scattering. 
Nevertheless, as represented by figure 3 in red, still 98 % of the electrons are contained within 
~ 5° of the vertical direction when they have reached the photosphere. Therefore, the vertical 
injection employed here is valid as long as the magnetic fields are sufficiently normal to the solar 
photosphere, the electrons' pitch angle is initially concentrated near ~ 0°, and non-Coulombic 
(e.g., magnetohydro-dynamic, or due to magnetic mirroring) pitch-angle scatterings can be 
neglected. 

The bremsstrahlung photons emitted by the almost vertically moving electrons are 
strongly forward collimated, as indicated by figure 7. Therefore, most of them are expected 
to enter the atmosphere and undergo Compton scattering or photoabsorption. Although some 
of these photons will die, others will be scattered back to escape out of the photosphere. We 
collected these outcoming photons at different viewing angles f3 (figure 1), and averaged the 
results over < ip < 2ir. Figure 11 shows the photon spectra obtained in this way for equal 
intervals in cos/3 (hence over an equal solid angle). 

In figure 11, spectra with small (3, simulating disk-center flares, exhibit a clear "knee", 
which represents the effect of single Compton scatterings. In the case of (3 ~ (blue in figure 
11), the knee energy becomes ~ 1/2 m e c 2 , because of the strong energy degradation in the 
single Compton back-scatterings. Photons with energies above the knee result either directly 
from the bremsstrahlung source, or through multiple Compton scatterings with a relatively 
small energy loss in each step. The spectrum in this region has a much steeper slope than in 
the region below the knee. 

As j3 gets larger in figure 11, the knee energy increases, because the angle needed in the 
single Compton scattering decreases. Nevertheless, the spectrum remains much softer (e.g., 
7 ~ 3 in the 0.2-0.6 MeV range) than is expected when the Compton effects are not taken into 
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account (7 ~ 1.3; figure 6b). The observed photon flux below the knee energy decreases as (3 
increases and approach 90°, due to the same limb darkening effect as in visible light. 

3.3. Slant injection 

As we have confirmed in figure 11, bremsstrahlung photons back-scattered from the solar 
photosphere exhibit very steep spectra as a result of the Compton degradation. In order to 
explain the origin of the observed flare spectra often extending beyond 1 MeV, we must then 
considered photons which reach us after Compton scatterings with much smaller angles. For 
this purpose, we have to consider electrons penetrating into the atmosphere with shallow angles 
(i.e., large a), namely "slant injection". 

We then examined the case of slant injection, by shooting 10 6 power-law electrons (in the 
1-100 MeV energy range), into our realistic solar model. Specifically, we repeated the Geant4 
simulation under the same condition, but assuming this time an injection angle of a = 80° 
(figure 1), instead of the vertical injection (a = 0°) considered in the previous subsection. We 
now inject the electrons from a height of 100 km (i.e., neglecting the uppermost two zones in 
figure 10), to avoid them escaping sideways before arriving at thicker regions. Figure 12 shows 
the results of this simulation for the same three intervals of (3 as figure 11, again averaged 
over the observer's azimuth angle </?. The spectrum becomes considerably harder than in the 
case of vertical electron injection, and extends much beyond ~ 1 MeV when (3 is rather large 
(0.1 < cos/3 < 0.2). In this case, electrons make grazing incidence on the atmosphere, and some 
of their bremsstrahlung photons get nearly forward scattered to leave the zone with a small 
energy loss, thus producing the hard power-law continuum. 

In this way, we have Monte-Carlo simulated slant injections for a family of electron 
spectra with various indicies 5 (again distributed over 1-100 MeV), assuming two representative 
injection angles of a = 60° and a = 80°. Figure 13 summarizes the photon index 7 of 1-10 MeV 
gamma-rays, collected over the most favorable viewing angle of 0.1 < cos/? < 0.2 (i.e., rather 
grazing to the solar photosphere) and integrated over the observer's azimuth. Thus, even with 
the extreme assumption of a = 80°, the simulated gamma-ray spectra become significantly 
softer (typically by ~ 1 in 7) than the original bremsstrahlung spectra before the photons are 
Compton scattered (reproduced in figure 13 in blue). Moreover, 7 further increases by ~ 0.5 
as a decreases from 80° to 60°. These results imply that the Compton scattering significantly 
soften the emergent photon spectra, of which the effect depends sensitively on a and j3. 

An interesting prediction of figure 13 is that the observed gamma-ray photon index 
should appear in a relatively narrow range of 7 = 1.7 — 2.6, regardless of the electron index 6. 
Of course, we would observe considerably larger values of 7 when either a or (3 is rather small. 
However, such a spectral component would not be easily detected in the MeV energy region. 

Although our simulations of the slant injection neglected the uppermost two coronal 
zones, their effects may not be necessarily negligible, because the large injection angle a = 80° 
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(or 60°) effectively increases the electron path length in each zone by (cos80°) _1 = 5.8 [or 
(COS60 )" 1 = 2.0] times. Since the electron energy losses in these two regions, if properly taken 
into account, would reduce the maximum electron energies while flatten the low-energy electron 
slope (figure 3), the emergent photon flux above ~ 1 MeV would further be suppressed if the 
electron injection in actual flares takes place at a height much exceeding the assumed 100 km. 
In other word, our simplification employed in simulating the slant injection makes the case 
more conservative. 
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4. Discussion 



4-1. Summary of simulations 

In order to investigate the effects of Compton scatterings on the solar gamma-ray spec- 
tra, we have conducted Monte Carlo simulations of energetic electrons, impinging on the solar 
atmosphere, using Geant4 as our basic toolkit. We modeled the solar atmosphere as vertically 
stacked parallel hydrogen zones, each assumed to be uniform. The density gradient was ex- 
pressed by assigning higher densities to lower zones. We injected electrons which are power-law 
distributed over the 1-100 MeV range, changing their spectral index 5 and incident angle a. 
The results of the present simulations can be summarized as follows. 

1. The observed photons hardly reach MeV energies if the electrons are injected vertically, 
because the bremsstrahlung-produced gamma-rays must be Compton back-scattered to 

reach the observer, and hence their energies are subject to strong Compton degradation. 

2. The emergent gamma-ray spectrum, integrated over the observer's azimuth, gradually 

hardens as the injection angle a increases (becoming more "slant"), and also as the ob- 
server's polar angle (3 gets larger. 

3. As long as a and f3 are both relatively large, the azimuthally-averaged 1-10 MeV gamma- 
ray spectra exhibit 7 = 1.7 ~ 2.6 over a relatively wide range of 5. 

These results generally agree with, and further extend, the previous works by, e.g., Bai & 
Ramaty (1978) and Kontar et al. (2006). 

4-2. Effects ignored in the simulation 

There are some effects which are not considered in our simulation. Among them, those 
of the secondary electrons and helium were already described in section 2.2 and section 3.1, 
respectively. As the electron energy increases, their synchrotron energy loss becomes significant. 
However, in the chromospheric region in which we are interested, the magnetic field intensity 
is typically ~ 100 G, and hence the rate of synchrotron loss of 100 MeV electrons is only ~ 1 
% of their Coulomb loss, and only ~ 5 % of their bremsstrahlung loss. This is because the 
Coulomb and bremsstrahlung losses are rather large due to the high matter density. Thus, the 
synchrotron process can be neglected. 

We also ignored the inverse Compton process between flare-accelerated electrons and 
ambient photons, because electrons with the maximum energy assumed, 100 MeV, can produce 
at most hard X-rays up to ~ 100 keV by scattering off the solar visible photons. However, if 
the electrons have a very flat spectrum well extending to energies above 350 MeV, the inverse 
Compton process may become a key process in creating gamma-ray photons. 

The size of the world used in our simulation, 10000 km x 10000 km x 20000 km, 
simulates a typical size of a foot point of a magnetic tube. The number of photons which 
escaped from the boundaries below the photosphere is less than 10 -4 of that of produced 
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photons, and energies of these escaped photons are less than 100 keV each. Hence we ignored 
those photons in the simulation. Possible effects of the relatively low heights of our slant 
injection was already considered in subsection 3.3. 

In the present simulation, we truncated the power-law electron spectrum at 100 MeV. 
In order to examine whether this abrupt cutoff produces any artifact, we repeated some of the 
simulations by changing the maximum cutoff energy to 1 GeV. q Then we have found that 7 
decreases only slightly, typically 0.1 or less, if 8 is larger than 2. 

Finally, a separate gamma-ray component could be produced without the Compoton 
effect, simply if the electron spectrum consists of two power-law components with different 
slopes. Even so, the good agreement between the predicted range of power-law indices and the 
Yohkoh data suggests that the Compton scattering is playing a significant role (see the next 
section). 

4-3. Comparison with observations 

In order to evaluate the significance of the results of our simulation, let us compare them 
with actual observations of hard X-ray and gamma-ray emission from solar flares. For that 
purpose, we have also to relate the simulation parameters a and /3 to observational parameters, 
including pitch-angle distributions of electrons, the tilt of magnetic field lines, and the flare 
longitude on the solar disk. As the observational data, we utilize 40 flares that show significant 
gamma-ray emissions, selected out of 2788 X-ray flares observed from 1991 October to 2001 
December by Yohkoh (Matsumoto et al. 2005). Yohkoh, the Japanese solar satellite launched in 
1991, was able to take four-color hard X-ray flare images with the Hard X-ray Telescope (HXT) 
(Kosugi et al. 1991; Kosugi et al. 1992) and perform spectro-photometry over the 0.2 — 30 MeV 
range with the Gamma-Ray Spectrometer (GRS) (Yoshimori et al. 1991). 
4-3.1. Vertical magnetic fields 

As the simplest case, we may consider a condition where the magnetic field lines are 
perpendicular to the photosphere. In this case, we can identify a with the pitch-angle of an 
electron, and the effect of electron gyration around the magnetic field lines is fully represented 
by our photon collection method which takes an average over ip. If the electrons have a narrow 
pitch-angle distribution around 0, the strong Compton degradation, described in section 3.2 
and shown in figure 11, will make the observed photon spectrum very steep, and will prevent 
it from reaching MeV energies. If, on the other hand, the pitch-angle distribution is rather 
broad, electrons with larger values of a will emit harder spectra, as represented by figure 12. 
Conversely, as long as the field lines are close to vertical, flares with significant gamma-ray 
emission must involve broad pitch-angle distributions of electrons. 

Analyzing the 40 gamma-ray flares detected by Yohkoh, Matsumoto et al. (2005) found 
that their gamma-ray (typically 1 MeV) to hard X-ray (~ 70 keV) flux ratios scatter largely (by 
an order of magnitude), with no correlations to their hard X-ray spectral slopes. Hence they 
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conclude that there is a hidden parameter which causes the gamma-ray flux to vary significantly 
relative to the hard X-ray extrapolation, in such a way that the gamma-ray component often 
appears as a separate spectral hard tail (see their figure 6). Given the present results, this 
parameter is very likely to be the electron pitch-angle distribution, with broader distributions 
yielding higher gamma-ray to hard X-ray flux ratios. 

In addition to the electron injection angle a, the viewing angle (3 has been confirmed 
to strongly affect the flare gamma-ray spectra (section 3). Assuming again the vertical field 
geometry, we can identify f3 with the angle 9 V which our line-of-sight to the flare makes relative 
to the normal to the photosphere. Hereafter, we simply call 9 V "flare- viewing angle". Using 
the flare longitude / and latitude b, it can be calculated as 



To examine the effect of (3, we plotted in figure 14 the gamma-ray to hard X-ray hardness 
ratio of the 40 Yohkoh flares, as a function of 8 V calculated via equation (12). It is similar to 
figure 7 of Matsumoto et al. (2005), but differs from it in that 9 V is used instead of \l\. As 
predicted by our simulation, the hardness ratio clearly increases with 9 V ; this effect was already 
known as "gamma-ray limb brightening" in solar flares using SMM data (Vestrand et al. 1987). 
Vestrand et al. (1987) interpreted this effect in terms of the radiation anisotropy: supposing 
that the electrons have a pitch-angle distribution which increases with angle from the outward 
normal, then higher-energy continua are expected to be more strongly beamed parallel to 
the photosphere, because harder bremsstrahlung photons are more strongly forward-beamed 
than softer ones even when emitted by electrons of the same energy. We presume that this 
mechanism can explain at least part of the observed limb-hardening effects, while the Compton 
back-scattering process must be enhancing them. In particular, our scenario can explain limb 
hardening even when the electrons are mostly directed downwards. 

For a more quantitative comparison of our simulation with the Yohkoh observations, 
let us limit ourselves to limb flares which have absolute longitudes larger than 60°, because 
gamma-ray data of disk flares are of poor photon statistics than those of limb flares as is clear 
from figure 14, and because these disk flares are more contaminated by nuclear gamma-ray 
lines (Matsumoto 2002) which are unrelated to the accelerated electrons. Figure 15 shows a 
scatter plot between two hardness ratios of the 40 sample flares, one calculated in the hard 
X-ray range (57 — 93 keV vs 33 — 57 keV) while the other in the gamma-ray range (1.43 — 6.21 
MeV vs 0.22 — 1.43 MeV). In this "color-color" plot, the dashed line represents the condition 
that the two colors imply the same photon index 7, of which values are given in the figure. 
The hard X-ray (or gamma-ray) photon index of a flare can be found by drawing a vertical (or 
horizontal) line from the data point, and read its intersection with the dashed line. 

From figure 15, we can see that the gamma-ray spectral index of the 40 flares is dis- 
tributed in a typical range of 1.7-2.5. This agrees very well with the prediction by our Monte- 




(12) 
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Carlo simulations (figure 13; subsection 3.3), on condition that the electrons in these flares 
have broad pitch-angle distributions beyond a ~ 60°. Further, in figure 15, a majority of flares 
are distributed to the left side of the single power-law locus. Such a concave-shaped wide-band 
spectrum can be naturally explained as a consequence of a broad pitch-angle distribution of 
electrons, because electrons with small pitch angles are expected to emit steep-spectrum pho- 
tons, while those with large pitch angles emit much harder spectra. A superposition of these 
spectral components with different slopes will result naturally in a concave spectrum just as is 
observed. Thus, as long as the magnetic lines are assumed perpendicular to the photosphere 
and the flare viewing angle 9 V is large, flares with significant gamma-ray signals are inferred to 
have rather broad pitch-angle distributions. 
4-3.2. Slant magnetic fields 

As a more general case, we may briefly consider tilted magnetic-field conditions. In 
this case, it is somewhat difficult to directly compare the results of our simulation with the 
observational data, because the electron gyration can no longer be faithfully emulated by our 
method of taking an average over ip. Nevertheless, we may employ the present results to predict 
some aspects of that case in a qualitative way. 

Let us first consider limb flares as an extreme case. In this field-line tilt on the 

sky plane is considered to have little effects, because the most efficient gamma-ray emission (as 
seen from the observer) occurs over the gyration phase wherein the electrons are approaching 
us, and this phase is the same even the filed lines are tilted in that way. In contrast, a field-line 
tilt along our line-of-site is expected to have a prominent effect: if the field lines are tilted 
away from us, the observed spectra tend to reach higher energies, e.g. a few MeV, even if 
the electrons have a rather narrow pitch-angle distribution. If the lines are tilted toward us, 
the gamma-ray spectra are expected to soften, on the contrary. Then, the observed scatter in 
the gamma-ray to hard X-ray intensity ratio for limb flares (figure 14, with 9 V > 60°) may be 
contributed by the field-line tilts, in addition to the pitch-angle distribution. 

The other extreme case is a disk-centered flare, where the field-line tilt is specified only 
by its angle 9t with respect to the photosphere normal (i.e., our line-of-sight) without depending 
on the tilt direction. Then, an electron with a fixed pitch angle, cto is expected to sweep a range 
of injection angle, from a_ = |a — &t I to a + = «o + ^- Therefore, a larger value of 9 t is favorable 
in delivering strong gamma-ray signals to the observer, even though j3 is still limited to the 
unfavorable value of ~ 90°. If «o and 9 t are both rather large so that ct + > 90°, the electron is 
implied to be approaching us over a certain phase of its gyration, thus allowing a considerable 
fraction of the bremsstrahlung photons to directly reach us without being Compton scattered. 
Therefore, a disk flare may become a strong gamma-ray source under a limited condition, that 
the electrons have a broad pitch angle and the magnetic field lines are significantly tilted. 
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4-3.3. Case of a broad pitch- angle distribution 

Now that broad pitch-angle distributions of accelerated electrons are considered to be an 

important element of observing intense gamma-ray emission, it may be imperative to simulate 

such a case, and compare the results with gamma-ray spectra from actual solar flares. Therefore, 

we repeated our solar simulation, assuming that the electrons have a continuous incident-angle 

distribution as 
3 

oc — (1 — cos 2 a) , (13) 

which corresponds to M2(cosa) in Dermer & Ramaty (1986) except for covering only the 
< cos a < 1 range. When the magnetic fields are vertical, equation (13) becomes identical to 
the pitch-angle distribution of the accelerated electrons. Under this angular distribution, we 
injected 10 7 electrons in the same way as in section 3, assuming the incident electron spectrum 
to have a canonical index of 5 = 2.0, and to extend again over 1-100 MeV. Collecting the 
produced gamma-rays over a viewing angle range of 0.1 < cos/3 < 0.2, and convolving the 
results with an approximate response of the Yohkoh GRS (Yoshimori et al. 1991), we have 
obtained the spectrum shown in figure 16 in red. 

To compare with this simulation, we have chosen a GOES X3.7 class flare of 1998 
November 22, 06h37m UT, which is one of the gamma-ray brightest events among the 40 
Yohkoh flares (Matsumoto et al. 2005). Since this is an extreme limb flare located at S32W90, 
our viewing angle (3 should be close to 90°, in agreement with the simulation condition. Actually, 
the HXT highest-band images of this flare reveal intense hard X-ray emission from a pair of 
magnetic loop foot points, both located right on the solar west limb (Matsumoto 2002). The 
GRS data, accumulated over a period of 06:37:22-06:40:22, are shown in figure 16 in blue, 
after subtracting background (accumulated for 06:56:14-06:59:02), but without removing the 
detector response. Thus, the spectrum is featureless, and the 2.23 MeV neutron capture line 
is undetectable with its equivalent width being < 10 keV (Matsumoto 2002). Therefore, this 
flare is considered to be a typical "electron dominated" event. 

In figure 16, the simulation and the actual data agree reasonably well, at least over the 1- 
10 MeV band. When the GRS response is considered, these spectra can be approximated both 
by T = 2.0 ±0.1 (Matsumoto 2002), which is typical in figure 15. Incidentally, we would obtain 
from figure 13 V ~ 1.7 for the assumed 5 = 2.0, if we were observing direct thick-target emission. 
We should instead find T = 2.2 (for a = 80°) to V = 2.5 (for a = 60°), if the Compton back- 
scattering dominates. Therefore, the values of T ~ 2.0 found in figure 16 are in between these 
two cases. Indeed, an inspection of the Monte-Carlo photons revealed that ~ 60% of them are 
Compton-scattered gamma-rays, while the remaining ~ 40% are direct bremsstrahlung events. 
Given the good agreement between the two data sets, it is possible that the actual gamma-rays 
from the November 22 flare are also a mixture of the direct and scattered photons. 



19 



5. Conclusion 

In order to quantitatively estimate the spectra of energetic electrons generated in solar 
flares based on the observed gamma-ray spectra, we numerically studied the electron transport 
and gamma-ray emission in the solar atmosphere. As elementary processes, we mainly con- 
sidered Coulomb scattering, bremsstrahlung and Compton scattering. We modeled the solar 
atmosphere with a vertical stack of parallel zones with different densities, and neglected the 
magnetic fields. 

As a result of our simulation, we have found that the thick-target gamma-ray spec- 
tra emitted from magnetic loop foot-point can hardly reach MeV energies, because of heavy 
Compton degradation, unless the electrons make rather grazing angles to the photosphere and 
we are observing forward-scattered gamma-rays. When this Compton effect is taken into ac- 
count, the traditional 5 vs. 5 relation is drastically modified. The pitch- angle distribution of 
flare-accelerated electrons are suggested to be broad and extend to large angles for those flares 
of which the emitted photon spectra reach MeV energies. Actually, a gamma-ray spectrum 
simulated under electron injection with a broad pitch angle distribution reproduced reasonably 
well the Yohkoh spectrum of the 1998 November 22 flare. 

In terms of actual flares, the above results can be translated into the following conclu- 
sions. 

1. When the magnetic field lines are perpendicular to the photosphere, and the electrons 
have a narrow pitch- angle distribution around 0, the gamma-ray spectra are very soft and 
cannot reach MeV energies. 

2. Even in the same field configuration, the gamma-ray spectra become harder if the electrons 
have larger pitch angles. Therefore, the electrons producing gamma-rays are inferred to 
have relatively large pitch angles. 

3. Limb flares are expected to have harder spectra than disk flares, because of their smaller 
angles of Compton scattering. This agrees with the observed fact that the gamma-ray to 
hard X-ray intensity ratio increases as the flare longitude gets closer to the solar limb. 

4. In case of a broad pitch angle distribution, electrons with small pitch angles emit steep 
hard X-rays, while those with large angles emit photons with flatter spectra. This can 
explain another Yohkoh result, that the gamma-ray slope tends to be flatter than the 
contemporaneous hard X-ray slope. 

5. Under the conditions of 5 = 1 — 3 and a > 60°, the values of 7 obtained in our simulations 
in the 1-10 MeV range are concentrated over 1.7-2.5. This agrees again with the Yohkoh 
observations. 

6. A disk flare may become a strong gamma-ray source under a limited condition, that the 
electrons have a broad pitch-angle and the magnetic field lines are significantly tilted. 

7. The gamma-ray spectra of actual solar flares are expected to be an appropriate mixture 
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of the direct and Compton-scattered photons. 
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Fig. 1. A coordinate system in the present simulation. The normal to the solar photo- 
sphere is z axis, and the electron injection angle measured from +z direction is denoted 
a. The viewing direction is specified by a polar angle (3 and an azimuth angle ip. 
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Electron Energy (MeV) Electron Projected Angle [rad] 



Fig. 2. A Geant4 simulation of electrons with an initial kinetic energy of 50 MeV, injected vertically into 
a hydrogen target with a uniform density of 3.2 x 10~ 7 g cm~ 3 (the simple solar atmosphere). Data are 
collected at depths of 100 km (red), 200 km (blue), and 300 km (green), (a) The energy spectrum, (b) The 
distribution (per unit projected angle) of angular deflection, measured from the initial direction of injection. 
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Fig. 3. The same as figure 2, but when the initial electrons are distributed between 1 
and 100 MeV with a power-law spectrum of photon index 5 = 1.2. Red. blue, green, 
cyan, and magenta indicate depths of 100, 200, 300, 400 and 500 km, respectively. 
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Fig. 4. Thin-target bremsstrahlung photon spectra, produced in a hydrogen gas slab (the thin 
atmosphere model) by mono-energetic electrons with kinetic energies of 10 MeV (magenta) and 
50 MeV (red). Crosses show the Geant4 simulation results, using the thin atmosphere 
model geometry, while solid curves refer to the analytic expression by Koch & Motz (1959) 
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Emitted Photon Energy (MeV) 

Fig. 5. Thin-target brcmsstrahlung spectra produced by incident electrons having a power-law index of 
8 = 0.8 and 1.5, distributed over the 1-100 MeV range. Crosses show the Geant4 simulation results using 
the thin atmosphere model geometry, while solid curves refer to the analytic expression by equation (6) 



26 



Incident powr-law (1-100 MeV) 



Incident powr-law (1 -1 00 MeV) 




1.5 2 
Incident Electron Index 




1.5 2 
Incident Electron Index 



Fig. 6. (left) Relation between the index 8 of power-law distributed (1-100 MeV) electrons and the 
resulting thin-target bremsstrahlung photon index 7. The red lines is the calculation from equation 
(6), blue crosses show the results of Geant4 simulations, and the black line refers to the non-rela- 
tivistic formula represented by equation (3). The photon index refer to the 1-10 MeV range, (right) 
The same as left figure, but for the thick target emission. The analytical results (red) is calcu- 
lated from equation (10), while the blue line represents the non-relativistic relation of in equation (2). 
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Fig. 7. Angular distributions (per unit polar angle measured form the electron motion) of thin- 
target bremsstrahlung photons, produced by electrons with the kinetic energies of 10 MeV. Crosses 
represent the Monte Carlo simulations, while solid curves are analytic expressions by equation (7). 
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Fig. 8. The same as figure 4, but obtained under the thick-target condition using the 
thick atmosphere model. The solid curves represents analytic calculations by equation (9). 
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Fig. 9. The same as figure 5, but obtained under the thick-target condition using the 
thick atmosphere model. The solid curves represents analytic calculations by equation (10). 
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Fig. 10. Density profiles of our "solar atmosphere" (blue) compared with the analytical approxima- 
tion by equation (11), in the regions above the photosphere (z > 0) panel (a) and below the photo- 
sphere (z < 0) panel (b). Green curves show the column density integrated downwards from z = +oc. 



31 




0.1 0.2 0.3 0.4 
Photon Energy (MeV) 

Fig. 11. Monte-Carlo simulated spectra of gamma-rays emergent from our "solar atmosphere", 
when 1-100 MeV electrons with S = 1.2 are vertically injected into it. The photons, 
mostly Compton back-scattered, are averaged over the observer's azimuth ip (figure 1), and 
collected at different viewing angle /3; red, green, and blue crosses represent the pho- 
tons satisfying 0.1 < cos (3 < 0.2, 0.5 < cos (3 < 0.6, and 0.9 < cos /? < 1.0, respectively. 
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Fig. 12. The same as figure 11, but simulating a highly slant injection with a = 80° instead of the 
vertical injection. Note that the displayed energy range is different from that used in figure 11 
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Fig. 13. Photon indicics 7 of the 1-10 McV gamma-rays, simulated under two favorable slant-injec- 
tion conditions, shown as a function of the electron power-law index 6 (over 1-100 MeV). Red aster- 
isks and green crosses represent a = 80° and a = 60°, respectively. In both cases, the emergent pho- 
tons are accumulated over 0.1 < cos/3 < 0.2 and < ip < 2ir. Blue crosses are the same as figure 6b, 
namely the spherically integrated bremsstrahlung photon index before subjected to the Compton process. 
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Fig. 14. The gamma-ray (0.3-13 McV) vs. hard X-ray (53-93 kcV) peak count 
tios of the 40 gamma-ray emitting solar flares observed by Yohkoh (Matsumoto et 
2005), plotted as a function of their flare viewing angles calculated via equation 



35 



0.5 



~l 1 1 1 1 1 1- 



I ong i tude | > 60 deg, 

r=i.5> 



u 

Oh 

^ n 9 



0.1 



/ r=2.5 



0.05. 
0. 



r=3.o 



_i i i i i i i L 



0.5 1 

H/M2 



Fig. 15. Color-color plots of the 40 gamma-ray emitting solar flares with absolute longitudes 
> 60° , observed with Yohkoh. Abscissa is the hard X-ray color between the H-band (57 — 93 
kcV) and M2-band (33 — 57 kcV) of the HXT, while ordinate is the gamma-ray color between 
the PC2-band (1.43—6.21 McV) vs. PCl-band (0.22—1.43 McV) of the GRS. The dashed 
line indicates the locus on which the two colors are represented by a single photon index 7. 
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Fig. 16. A gamma-ray spectrum (blue) of the 1998 November 22 solar flare, acquired with the 
Yohkoh GRS over a time period of 06:37-06:40 UT. It is compared with a simulation (red crosses), 
computed assuming that 10 7 electrons with a spectrum of 8 = 2.0 (1-100 MeV) and an inci- 
dent angle distribution of equation (13) are injected, and the emergent photons are observed with 
0.1 < cos (3 < 0.2. The simulated spectrum has been convolved with an approximate response of the 
GRS. The normalization of the simulated spectrum is arbitrarily adjusted to match the observed data. 
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